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Abstract. We use ideas from algebraic geometry and dynamical systems to explain some 
' ways that control points influence the shape of a Bezier curve or patch. In particular, we 

establish a generalization of Birch's Theorem and use it to deduce sufficient conditions on 
the control points for a patch to be injective. We also explain a way that the control points 
influence the shape via degenerations to regular control polytopes. The natural objects of 
this investigation are irrational patches, which are a generalization of Krasauskas's toric 
. patches, and include Bezier and tensor product patches as important special cases. 

< 

' Introduction 

The control points and weights of a Bezier curve, Bezier patch, or tensor-product patch 
CSJ ■ govern many aspects of the curve or surface. For example, they provide an intuitive means 

^ ! to control its shape. Through de Castlejau's algorithm, they enable the computation 

^ I of the curve or surface patch. Finer aspects of the patch, particularly continuity and 

CN ■ smoothness at the boundary of two patches are determined by the control points and 

weights. Global properties, such as the location of a patch in space due to the convex hull 
^ i property, also depend upon the control points. When the control points are in a particular 

QQ I convex position, then the patch is convex 

O ■ We apply methods from algebraic geometry, specifically toric geometry, to explain how 

>• . some further global properties of a patch are governed by the control points. We first 

I investigate the self-intersection, or injectivity of a patch. We give a simple and easy-to- 

^ ' verify condition on a set of control points which implies that the resulting patch has no 

self-intersection, for any choice of weights. For 3-dimensional patches as used for solid 
modeling, injectivity is equivalent to the patch properly parameterizing the given solid. 
This uses Craciun and Feinberg's injectivity theorem [1] from the theory of chemical reac- 
tion networks, which may be seen as a generalization of Birch's Theorem from algebraic 
statistics. 

A second global property that we investigate is how the shape of the patch is related 
to the shape of a control polytope. This is a piecewise linear triangulated surface whose 
vertices are the control points. It is regular if the underlying triangulation comes from a 
regular triangulation of the domain polytope of the patch. We show that regular control 
polytopes are the limits of patches as the weights undergo a toric deformation corresponding 
to the underlying regular triangulation, and that non-regular control polytopes can never 
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be such a limit. This gives a precise meaning to the notion that the shape of the control 
net governs the shape of the patch. 

This line of inquiry is pursued in terms of Krasauskas's toric patches [H], as it relies 
upon the structure of toric varieties from algebraic geometry. The correct level of gener- 
ality is however that of irrational (toric) patches, which are analytic subvarieties of the 
simplex (realized as a compactified positive orthant) that are parameterized by monomials 
where a; is a vector of positive numbers and the exponent vector a has real-number 
coordinates. (This is a usual toric variety when a has integer coordinates.) While irra- 
tional patches may seem exotic for modeling, they occur naturally in statistics as discrete 
exponential families [2] and their blending functions may be computed using iterative pro- 
portional fitting (IPF) [7], a popular numerical algorithm from statistics. Furthermore, 
these blending functions have linear precision. For toric patches, this was observed in [18] 
and developed in [10], and the analysis there carries over to irrational patches. While we 
work in this generality, our primary intent (and the main application) is to shed light on 
properties of Bezier curves, surfaces, and 3-dimensional patches. 

We recall the standard definition of a mapping via control points and blending functions, 
and then the definitions of toric Bezier patches in Sect. [U There, we also illustrate some of 
our results on examples of Bezier curves. In Sect. [21 we introduce irrational toric patches, 
recalling the geometric formulation of a toric patch and the use of iterative proportional 
fitting to compute these patches, explaining how these notions from [10] for toric patches 
extend to irrational patches. The next two sections contain our main results. We study 
injectivity of patches in Sect. [3], and discuss degenerations to control polytopes in Sect. HI 
Appendices A and B contain technical proofs of some theorems. 



We interpret the standard definition of a mapping via control points and blending func- 
tions (see for example [T2l §2]) in a general form convenient for our discussion. All functions 
here are smooth (C°°) where defined and real- valued. Let M> be the set of strictly positive 
real numbers and ]R> the set of non-negative real numbers. We will use the following typo- 
graphic conventions throughout. Vector constants (control points, indexing exponents, and 
standard basis vectors) will be typeset in boldface, while vector variables will be typeset 
in standard math italics. 

Let ^ be a finite set of points that affinely span M'^, which we shall use as geometrically 
meaningful indices. A control point scheme for parametric patches, or (parametric) patch, 
is a collection /5 = {/3a | a G A} of non- negative functions, called blending functions. The 
common domain of the blending functions is the convex hull A of A, which we call the 
domain polytope. We also assume that the blending functions do not vanish simultaneously 
at any point of A, so that there are no basepoints. 

Lists B := {ba | a G A} C of control points and positive weights w := {w^ G M> | 
a G .4} G M> together give a map F: A — M" defined by 



1. ToRic Bezier Patches 





The denominator in f 1 1.1 1) is positive on A and so the map F is well-defined. 
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Remark 1.2. We will refer to both Pa{x) and WaPa.{x) as the blending functions of a 
patch. This generality of separating the weights from the blending functions will be used 
in Sect. Hlwhen we investigate the effect of systematically varying the weights of a patch 
while keeping the control points and blending functions constant. 

The control points and weights affect the shape of the patch which is the image of the 
map F (II. ip . For example, the convex hull property asserts that the image F{A) of the 
patch lies in the convex hull of the control points. To see this, note that if we set 

then l3g^{x) > and 1 = X]ae^'^a(^)- Then formula (11. ip becomes 

Fix) = , 

SO that F{x) is a convex combination of the control points and therefore lies in their convex 
hull. In fact, if there is a point x G A at which no blending function vanishes, then any 
point in the interior of the convex hull of the control points is the image F{x) of some 
patch for some choice of weights. In this way, the convex hull property is the strongest 
general statement that can be made about the location of a patch. 

Another well-known manifestation of control points is the relation of a Bezier curve to 
its control polygon. Fix a positive integer m and let ^ := | i = 0, . . . , m} so that A is 
the unit interval. The blending functions of a Bezier curve are the Bernstein polynomials, 

(3^{x) (= /3^{x)) := (7)x*(l-x)™-^ . 

The control polygon of a Bezier curve with control points bo, bi, . . . , is the union of the 
line segments bo,bi, bi,b2, bm-i,bm between consecutive control points. Figured] 
displays two quintic plane Bezier curves with their control polygons (solid lines). The 




Figure 1. Quintic Bezier curves. 

convex hulls of the control points are indicated by the dashed lines. The first curve has no 
points of self-intersection, while the second curve has one point of self-intersection. While 
this self-intersection may be removed by varying the weights attached to the control points. 
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by Theorem 13.71 it is impossible to find weights so that a curve with the first set of control 
points has a point of self-intersection. 

We will also show that the control polygon may be approximated by a Bezier curve. We 
state a simplified version of Theorem 14.41 from Sect. HI 

Theorem. Given control points in M" for a Bezier curve and some number e > 0, there 
is a choice of weights so that the image F[0, 1] of the Bezier curve lies within a distance e 
of the control polygon. 

In Figure [2], we display one of the quintic curves from Figure [H but with weights on 
bo— bs of (1, 20^, 20^, 20^, 20^, 1) and (1, 300^, 300^, 300^, 300^, 1), respectively. The control 
polygon for the second curve is omitted, as it would obscure the curve. The first curve 
lies within a distance e = 0.13 of the control polygon and the second within a distance 
e = 0.02, if the control polygon has height 1. 




Figure 2. Degenerating quintics. 

1.1. Toric Patches. Krasauskas [14j introduced toric patches as a generalization of the 
classical Bezier and tensor product patches. These are based upon toric varieties from 
algebraic geometry and their shape may be any polytope with integer vertices. The arti- 
cles [HI [18] provide an introduction to toric varieties for geometric modeling. 
A polytope A is defined by its facet inequalities 

A = {x G M"' I < hi{x) ,i = !,...,£} . 

Here, A has i facets (faces of maximal dimension) and for each i = 1, . . . hi{x) = Vj-x+Q 
is the linear function defining the ith facet, where Vj G is the (inward oriented) primitive 
vector normal to the facet and q G Z. 

For example, if our polytope is the triangle with vertices (0, 0), (m, 0), and (0, m), 

(1.3) ?Ti^ := {{x,y) G R^\0<x,y, and < m - {x + y)} , 

then we have hi = x, h2 = y, and hs = m — x — y. Here, m ^ is the unit triangle ^ with 
vertices (0, 0), (1, 0), and (0, 1) scaled by a factor of m. 
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Let A G AnZ'^ he any subset of the integer points of A which includes its vertices. For 
every a G .4, Krasauskas defined the toric Bezier function 

(1.4) p^{x) := /ii(x)^^(")/i2(x)^^(")---/i^(a;)^^(") , 

which is non- negative on A, and the collection of all /3a has no common zeroes on A. These 
are blending functions for the toric patch of shape A. If we choose weights w G M"^ and 
multiply the formula (11.41) by Wa, we obtain blending functions for the toric patch of shape 
{A,w). 

Example 1.5 (Bezier triangles). When A is a scaled triangle or a product of such tri- 
angles, (11.41) gives the blending functions of the Bezier patch or Bezier simploid |B] with 
the corresponding shape. To see this for the scaled triangle m\S. (II. 3p . note that given an 
integer point a = G m^, and weight the multinomial coefficient := , 

then the corresponding blending function is 

= .,!(r-.-,)! ^v(m-x-y)— . 

This is almost the bivariate Bernstein polynomial, which is obtained by substituting mx 
and my for x and y, respectively, and dividing by m™. (This has the effect of changing the 
domain from to the unit triangle ^.) 



2. Irrational Patches 

Krasauskas's definition ( 11. 4p of toric Bezier functions still makes sense if we relax the 
requirement that the points ^ C M*^ have integer coordinates. This leads to the notion 
of an irrational patch (as its blending functions are no longer rational functions), which 
provides the level of generality appropriate for our investigation. 

Let ^ C M"' be a finite collection of points and set A C M'^ to be the convex hull of A, 
which we assume is a full-dimensional polytope. We may also realize A as an intersection 
of half-spaces through its facet inequalities, 

(2.1) A = {x G M'^ I hi{x) > for each i = 1, . . . J} . 

Here, A has i facets with the ith facet supported by the affine hyperplane hi{x) = where 
hi{x) = \i ■ X + Ci with Cj G M and Vj an inward pointing normal vector to the ith facet 
of A. There is no canonical choice for these data; multiplying a pair (vj, q) by a positive 
scalar gives another pair defining the same half-space. 

Following Krasauskas, we provisionally define (irrational) toric Bezier functions {Pa- A — ; 
]R> I a G A} by the same formula as (11.41) . 

P.{x) := /ii(x)^^(")/i2(x)''^(")---/i^(x)''^(") . 

These are blending functions for the irrational toric patch of shape A. While these functions 
do depend upon the choice of data (vj, q) for the facet inequalities defining A, we will see 
that the image F{A) of such a patch given by weights and control points is independent 
of these choices. 
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2.1. Geometric Formulation of a Patch. We follow Sect. 2.2 of [lOj, but drop the 
requirement that our objects are algebraic. Let ^ C M'^ be a finite subset indexing a 
collection of blending functions {/9a: A ]R> | a G A}, where A is the convex hull of A. 
Let M"^ be a real vector space with basis {ca | a G .4}. Set M> C M"^ to be the points with 



non-negative coordinates and let be those points with strictly positive coordinates. 
For z = {z^\ae A) e M^, set ^ 2 := Eae^t ^a- The A-simplex, /S^ C M^, is the set 

ZS7^ := {2 G : = 1} • 
We introduce homogeneous coordinates for If 2; G M> \ {0}, then we set 

[^a I a G A\ := ^^{z^ I a G ^) G . 

The blending functions {/3a: A M> | a G A\ give a C°° map, 

/? : A — > 

X I — > I a G ^] . 

The reason for this definition is that a mapping F: A ^ MJ^ f lLlI) given by the blending 
functions /?, weights w, and control points B factors through the map /3 : A ^ To see 
this, first note that the weights w G ]R> act on /S/'^: If z = [z^ \ sl E A] E then 

(2.2) w.z := [wa^a I a G ^] . 
The control points B define the map vrg : M" via 

71 B : z = (^a I a G ^) I — > ^ 2;aba . 

Then the mapping F (11. ip is simply the composition 

(2.3) A ^ ^ ^ . 

In this way, we see that the image /5(A) C of A under the map P determines the 
shape of the patch -F(A) (11.11) . Internal structures of the patch, such as the mapping of 
texture, are determined by how j3 maps A to /3(A). For example, precomposing /5 with 
any homeomorphism of A gives blending functions with the same image in ^-^j but with 
a different internal structure. 

For an irrational toric patch of shape A with blending functions (11.41) . the image /?(A) C 
is independent of the choice of normal vectors. For this, we first define the map 
(^^: by 

(2.4) ■ (xi,...,Xd) I — ^ [x^ : aG^] . 

Let X_4 be the closure of the image of the map ipj^^. When A C Z'^, this is the positive 
part [HI §4] of the toric variety parameterized by the monomials of A. When A is not 
integral, we call X_4 the (irrational) toric variety parameterized by monomials in A. 
In Appendix |B] we prove the following theorem. 

Theorem 2.5. Suppose that A dR'^ is a finite set of points with convex hull A. Let (3 = 
{/3a I a G A} be a collection of irrational toric Bezier functions for A. Then /3(A) = X_a, 
the closure of the image of (f^. 
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We prove this by showing that the restriction of the map /3 to the interior A° of A factors 
through the map y^^^. 

By Theorem I2.5[ the image of the irrational toric blending functions for A depends 
upon A and not upon the choice of toric blending functions for A. Thus the shape of 
the corresponding patch -F(A) (11.11) depends only upon A, the weights w, and the control 
points B. However, the actual parameterization of by A, and hence of F{A) does 
depend upon the choice of toric blending functions for A. 

To ensure that the patch F{A) has shape reflecting that of A, we require that the 
map /3: A — i> be injective. This also guarantees that the patch F{A) is typically 
an immersion. In the context of irrational toric patches, this injectivity is guaranteed by 
Birch's Theorem from algebraic statistics. For a standard reference, see [U p. 168]. When 
A C Z"^, Birch's Theorem follows from general results on the moment map in symplectic 
geometry [9| §4.2 and Notes to Chapter 4, p. 140]. 

Theorem 2.6 (Birch's Theorem). Suppose A <Z Mf^ is finite and let j3 he a collection of 
toric Bezier functions. If we choose control points B to be the corresponding points of A, 
{ba = a I a G A}, then the composition 

is a homeomorphism onto A. 

By Birch's Theorem and Theorem l2.5[ any two sets /5, (3' of toric Bezier functions of shape 
A differ only by a homeomorphism h: A A of the polytope A, so that (3' = (3 o k. 
In fact h restricts to a homeomorphism on all faces of A. As we are concerned with the 
shape of a patch and not its internal structure, we follow Krasauskas' lead and make the 
following definition. 

Definition 2.7. A(n irrational) toric patch of shape A is any set of blending functions 
/? := {/3a : A — > M> | a G A} such that the map : A ^ X_4 is a homeomorphism. 

The projection map txb '■ IR*^ appearing in Birch's Theorem induced by the choice 

B = {ba = a I a G A} of control points is called the tautological projection and written tt^. 
Restricting the tautological projection to X_4 gives the algebraic moment map fi: X_4 —>■ A. 
The components of its inverse fi~^ : A ^—>- provide a preferred set of blending functions 
for the patch. When A C Z'^, these were studied in [l0l|T8], where they were shown to have 
linear precision, and that they may be computed by iterative proportional fitting (IPF), 
a numerical algorithm from statistics [7j. These same arguments apply to irrational toric 
patches — the preferred blending functions have linear precision and are computed by IPF. 

Any patch has unique blending functions with linear precision [TUl Theorem 1.11]. While 
the classification of toric patches for which these preferred blending functions are rational 
functions remains open in general, it has been settled for surface patches {d = 2) [TT], and 
this places very strong restrictions on higher-dimensional patches. 

2.2. Iterative Proportional Fitting for Toric Patches. In algebraic statistics, z^SZ"^ 
is identified with the probability simplex parameterizing probability distributions on data 
indexed by A. The image X^^^ := w.X^ of under the map (12.41) and translation 
by w ( 12. 2p is known as a toric model [Ml §1.2]. It is more common to call this a log-linear 
model, as the logarithms of the coordinates of ipj^^ are linear functions in the logarithms 
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of the coordinates of M^, or a discrete exponential family as the coordinates of <y9_4 are 
exponentials in the logarithms of the coordinates of M^. 

The tautological map appears in statistics as follows. Given (observed) normalized data 
q G the problem of maximum likelihood estimation asks for a probability distribution 
in the toric model, p G with the same sufficient statistics as q, nip) = ita{p) = T^Ail)- 

By Birch's Theorem, the point p G is unique and hence 

P = /^"^(7r.4(g)) ■ 

Thus inverting the tautological projection is necessary for maximum likelihood estimation. 

Darroch and Ratcliff [7] introduced the numerical algorithm of iterative proportional 
fitting, also known as generalized iterative scaling, for computing the inverse of the 
tautological projection. We now describe their algorithm. 

Observe first that the toric patch X_4 „, does not change if we translate all elements of A 
by a fixed vector b, (a i— > a + b), so we may assume that A lies in the positive orthant M>. 
Scaling the exponent vectors in ^ by a fixed positive scalar t G M> also does not change 
Xa,w as x (— > X* is a homeomorphism of M> which extends to a homeomorphism of M^. 
Thus we may assume that A lies in the standard simplex in M'^, 

Ad = {xGMi I Ea:< 1} . 

Lastly, we hft this to the probability simplex c M>"^\ 

:= {yGMfMEl/ = l} , 

by 

A 3 Si ^ a+ := a) G . 

Since for t G M> and x G M>, 

we see that replacing A by this homogeneous version also does not change X_4. 

We describe the algorithm of iterative proportional fitting, which is Theorem 1 in [7j. 

Proposition 2.8. Suppose that A C IS'^'^^ has convex hull A and q G /57-^. Set y : = 
7r_4(g) G A. Then the sequence of points 

{pM I m = 0,1,2...} C 

whose s.- coordinates are defined by p^^ := Wa and, for m > 0, 



(m+l) _ (m) . y 



converges to the unique point p G X_4 „, such that fi{p) = vr_4(p) = vr_4(g) = y. 

We remark that if A is not homogenized then to compute fi~^{y) for y G A, we first 
put A into homogeneous form using an affine map ip, and then use iterative proportional 
fitting to compute T^^+{^{y)) = T^A^iv)- We also call this modification of the algorithm of 
Proposition 12.81 iterative proportional fitting. Thus iterative proportional fitting computes 
the inverse image of the tautological projection. 
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3. Injectivity of Patches 

Birch's Theorem (Theorem 12.61) states that for one particular choice of control points, 
namely {bg := a | a G A} and all weights w G M^, the mapping F (11.11) of a toric 
patch of shape {A, w) is a homeomorphism onto its image. From this, we can infer that for 
most choices of control points and weights, this mapping is at least an immersion. To study 
dynamical systems arising from chemical reaction networks, Craciun and Feinberg prove 
an injectivity theorem for certain maps, which we adapt to generalize Birch's Theorem. 
This will give conditions on control points i3 C which guarantee that for any choice w of 
weights, the resulting mapping F (II. ip of a toric patch of shape {A, w) is a homeomorphism 
onto its image. This result has several consequences concerning the injectivity of toric 
patches. 

Let us first give the Craciun- Feinberg Theorem. Let Y = {yi, . . . , y^} C M" be a finite 
set of points which affinely spans M". For k G M™ and Z := {zi, . . . , z^} C M", consider 
the map : M!^ M" defined by 

m 

(3.1) Gk{x) := ^hx^^7.i . 

i=l 

Theorem 3.2 (Craciun- Feinberg). The map Gk is injective for every k G M> if and only 

if the determinant of the Jacobian matrix, 



Jac(Gfc) = 




does not vanish for any x G M" and any k G M> . 
We give a proof in Appendix \^ 

The condition of Theorem 13.21 that the Jacobian Jac(Gfc) does not vanish for any x G 
is reminiscent of the Jacobian conjecture [13], which is that the Jacobian of a polynomial 
map G: C" — s> does not vanish if and only if the map G is an isomorphism. Since 
we are restricted to x G M", it is closer to the real Jacobian conjecture, which is however 
false [T7], and therefore not necessarily relevant. 

The condition of Theorem 13.21 is conveniently restated in terms of Y and Z. For a list 
/ = {ii, . . . , in} C {1, . . . , m}, which we write as / G ('™') , let Yj be the determinant of the 
matrix whose columns are the vectors , . . . , y^^ , and define Zj similarly. In Appendix \^ 
we deduce the following corollary. 

Corollary 3.3. The map Gk (13. ip is injective for all k G M."^ if and only if {Yj ■ Zj) ■ (Yj ■ 

Zj) > for every I, J E ('^^) and at least one product Yj ■ Zj is non-zero. 

This leads to a generalization of Birch's Theorem. An ordered list pq, . . . ,pd of affinely 
independent points in M.'^ determines an orientation of M'^ — simply consider the basis 

Pi-Po, P2-P0, ■■■ , Pd-Po ■ 

Let A and B = {ba | a G A} be finite sets of points in M.'^. Suppose that {ao, . . . , a^} is an 
affinely independent subset of A. If the corresponding subset {bag, . . . , ba^} of B is also 
affinely independent, then each subset determines an orientation, and the two orientations 
are either the same or they are opposite. We say that A and B are compatible if either 
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every such pair of orientations is the same, or if every such pair of orientations is opposite. 
We further need that there is at least one affinely independent subset of A such that 
the corresponding subset of B is also affinely independent. Observe that compatibility is 
preserved by invertible affine transformations acting separately on A and B. 

In Fig. [3] shows three sets of labeled points. The first and second sets are compatible, 
but neither is compatible with the third. 



20 



30 



20 



10 



4« 



10 



40 



30 



10 



30 



20 



4 © 



Figure 3. Compatible and incompatible sets of points. 

We give our generalization of Birch's Theorem. Suppose that A C M"^ is the convex hull 
of A and {Pa- A M> | a G A} are toric Bezier functions for A. For any w E M>, let 
F^: A ^ M'^ be the toric patch of shape {A, w) given by the control points B C M'^: 

f3 4) F (x) - ^ag^^-^-(^)^- 

Theorem 3.5. The map F^ is injective for all w G M> if and only if A and B are 
compatible. 

As any set A is compatible with itself, this implies Birch's Theorem (Theorem 12. 6p . 

Example 3.6. Let l\ be the convex hull of {(0, 0), (1, 0), (0, 1)}. Set ^ := 3 l\ n and 
let w G Mt^ be the weights of a cubic Bezier patch (Example 11.51 with m = 3). We consider 
choices i3 C of control points that are compatible with A. For convenience, we will 
require that ba = a when a is a vertex and that if a lies on an edge of 3 then so does 
ba. For these edge control points, compatibility imposes the restriction that they appear 
along the edge in the same order as the corresponding exponents from A. The placement 
of the center control point is however constrained. We show two compatible choices of B 
in Fig. m On the left is the situation of Birch's Theorem, in which ba = a, and on the 
right we have moved the edge control points. The region in which we are free to move the 
center point is shaded in each picture. 

Theorem 1 5'. 51 Let (t, x) be coordinates for R'^^^ and consider the map : R^"^^ W^^^ 
defined by 

G^{t,x) = ^tx''Wa(l,ba) . 

aeA 

We claim that is injective if and only if is injective. 

Since F^ is the composition (12.31) A Xj^ w.X_a — ^ R'^, with the first map an 
isomorphism, F^ is injective if and only if the composition of the last two maps is injective. 
Since X_4 is compact, this will be injective if and only if its restriction to the interior 
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Figure 4. Compatible control points for the Bezier triangle. 

of is injective. Precomposing with the monomial parametrization (12.41) of X^, we see 
that Fw is injective if and only if the map '■ ffi> — > M"^ defined by 

is injective. 

Since G^it^x) = t ■ (7^^,(1, x), these values lie on a ray through the origin, and we invite 
the reader to check that this ray meets the hyperplane with first coordinate 1 at the point 
Thus is injective if and only if Hyj is injective, which is equivalent to the 
map Fw being injective. 

We deduce the theorem by showing that Gw is injective. This follows from Corollary 13.31 
as the condition that A and B are compatible is equivalent to 1/ ■ > for all / G {]^-^ , 

where F = {(1, a) \ae A}, Z = {(1, ba) | a G A} and we have m = ^ 
We now describe two applications of Theorem 13.51 to modeling. 

3.1. Solid Modeling with (Toric) Bezier Patches. In solid modeling, we represent a 3- 
dimensional solid by covering it with 3-dimensional patches, for example using Bezier toric 
patches as finite elements. Besides the obvious C° or higher continuity along the boundary 
as required, such Bezier finite elements should at least provide a one-to-one parametrization 
of their image, i.e. they should be injective. By Theorem 13.51 we may guarantee injectivity 
by requiring that the control points B be compatible with the exponents A. Moreover, if 
these sets are incompatible, then there is some choice of weights for which the patch is not 
injective. 

3.2. Injectivity of Bezier Curves and Surfaces. Typically, the exponents A and the 
control points do not lie in the same space; surfaces {A C M?) are modeled in 3-space 
{B C M^), or curves {A C M) in 2- or 3-space {B C or M^). Nevertheless, Theorem |33] 
gives conditions that imply injectivity of patches. 

Let p G M""*"^ be a point disjoint from a hyperplane, H. The projection M"'+^— H 
with center p is the map which associates a point x G M""''^ to the intersection of the line 
px with H. We use a broken arrow as the projection is not defined on the plane through p 
parallel to H. Identifying H with gives a projection map W^^^ > M". A coordinate 
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projection {xi, . . . ,Xn,Xn+i) ^ {xi, . . . ,Xn) is a projection with center at infinity. More 
generally, a projection M" > M*^ is a sequence of such projections from points. 

Theorem 3.7. Let A C M'^, w G M^, and B C be the exponents, weights, and control 

points of a toric patch and let A be the convex hull of A. If there is a projection ir : R" > 

M*^ such that A is compatible with the image 7t{B) of B, then the mapping F: A ^ 
given by the toric blending functions associated to A, the weights w, and control points B 
is injective. 

Proof. By Theorem 13.51 the composition vr o F is injective, from which it follows that F 
must have been injective. 

Example 3.8. For the curve on the left in Fig. [1] (which is reproduced below), the vertical 
projection maps the control points {bo, . . . , bs} to points on the line in the 

same order as the exponents A = {0,1,2,3,4,5}, which implies that the curve has no 
self- intersections. 




Figure 5. A compatible projection. 



4. Control Polytopes and Toric Degenerations 

The convex hull property asserts that the image, F{A), of a toric Bezier patch of shape 
{A, w) given by control points B = {ba | a G A} C M" and weights w G M> lies in the 
convex hull of the control points. When F{A) is a curve, the control points may be joined 
sequentially to form the control polygon, which is a piecewise linear representation of the 
curve. When F{A) is however a surface patch, there are many ways to interpolate the 
control points by triangles or other polygons to obtain a piecewise linear surface, called a 
control polytope, that represents the patch. The shape of this control polytope affects the 
shape of the patch. For example, when the control points have the form (a, A(a)) for A a 
convex function, then the patch is convex [6l[5]. Also, Leroy [15] uses a particular control 
polytope for the graph of a function to obtain certificates of positivity for polynomials. 

Among all control polytopes for a given set of control points, we identify the class of 
regular control polytopes, which come from regular triangulations of the exponents A. 
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These regular control polytopes are related to the shape of the patch in the following 
precise manner: There is a choice of weights so that a toric Bezier patch is arbitrarily close 
to a given control polytope if and only if that polytope is regular. 

4.1. Bezier Curves. It is instructive to begin with Bezier curves. A Bezier curve of degree 
m in M" with weights w is the composition (12. 3p . 

[0, 1] A, ^m+l ^ ^ 

where P = (Pq, . . . , Pm) with Pi{x) = (™')a;*(l — x)™""* for x E [0, 1]. Then the map P is given 
by Zi = - for i = 0, . . . , m. Here, (zq, • • • , ^m) G K>+^ with Zo + --- + Zrn = l 

are the coordinates for The image P[0, 1] C ^"^'^^ is defined by the binomials 

(4.1) (T)(7)^a^. - =0, for a + b = ^+J . 

To see this, suppose that (zq, . . . , Zm) G satisfies (14.11) . Setting x := zi/{mzo + zi), 

then we may solve these equations to obtain Zi = (™)x*(l — x)™"* = Pi{x). 
If w = {wo, . . . , Wm) G are weights, then w.p[0, 1] is defined in by 

(4.2) WiWj{"^){"l)zaZb - WaWb{'^){'^)ziZj = 0, for a + b = i + j . 
Suppose that we choose weights Wi := Dividing by f ("^-J)^ (14. 2 p becomes 

(4.3) (T)(7)^a^^ - t^^'"-''"-'\:){'^)z^z, = 0, for a + b = z+j . 

Since i + j = a + b, we may assume that a < i < j < b. Setting c := i — a = b — j > 1, we 
see that 

+ f -a^ -b'^ = c{i + a) - c{j + b) = c{i - j + a - b) < -2c < . 
If we consider the limit of these binomials (14. 3 p as t ^ oo, we obtain 

Za ■ Zb = if |a — 6| > 1 . 
These define the polygonal path in ^57™+^ whose ith segment is the edge 
{0^^^^^, X , 1 - X , 0^^^^^^ for a;G[0, 1], 

i—l m—i 

and whose projection to is the control polygon of the Bezier curve, which is the collection 
of line segments bo, bi, bi, b2, . . . , b^-i, b^. 

We illustrate this when m = 3 in Fig. El which shows three different Bezier curves having 
the same control points, but different weights Wi = for t = 1, 3, 9. In algebraic 

geometry, altering the weights in this manner is called a toric degeneration. The Bezier 
cubics are displayed together with the cubics w.p[0, 1] lying in the 3-simplex, which 
is drawn in M^. In these pictures, the projection vrg is simply the vertical projection 
forgetting the third coordinate. The progression indicated in Fig. [6], where the Bezier curve 
approaches the control polygon as the parameter t increases, is a general phenomenon. Let 
II • II be the usual Euclidean distance in M". 

Theorem 4.4. Suppose that Ft: [0, 1] ^ is a Bezier curve of degree m with control 
points B and weights Wi = fC™-*), ;= max{||ba|| : E B} . For any e > 0, if we 

have t > nm/e, then the distance between the control polygon and any point of the Bezier 
curve -Ff [0, 1] is less than e. 




Figure 6. Toric degenerations of a Bezier cubic. 

Proof. Let z E w.p[0,l] C where the weights are Wi = with t > nm/e. 

Suppose that h — a > 1 are integers in [0,m]. Then there exist integers i < j with 
a < i < j < b and a + b = i + j . By (14.31) . we have 

f'm\ (ra\ 

(T)(7) 

Since the binomial coefficients are log-concave, we have 



ZaZb — t ZiZj 



Using + — a? — b"^ < —2 and ;Zj + -Zj < 1, we see that 

1 

In particular, if |6 — a| > 1, then at most one of Za or z^ exceeds l/2t. We conclude that 
at most two, necessarily consecutive, coordinates of z may exceed l/2t. Suppose that i is 
an index such that Zj < l/2t if j ^ and let x := Zi^ihi^i + (1 — Zi_i)bj, a point 

along the zth segment of the control polygon. Since 

Em — 1 
> 1 ^ ' 

we have < 1 — Zi^i — Zi < ^^^^ < ^, and we see that 

Iki3(2;)-x|| = ^Zjbj - {zi^ihi^i + {1 - Zi_i)hi] 



< 



^illbjll + l^;^ - (1 - 2;i_i)|||bi 



m — 1 m Km 

< K + K — < = e. 

2t 2t t 
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This proves the theorem as 2; G w.p[0, 1] is an arbitrary point of the curve F[0, 1]. '^y' 

4.2. Regular Triangulations and Control Polytopes. In dimensions 2 and higher, 
the analog of Theorem 14.41 requires the notion of a regular triangulation from geometric 
combinatorics. Let ^ C M'^ be a finite set of points and consider a lifting function A : ^ ^ 
M. Let Pa be the convex hull of the lifted points 

A(^) := {(a, A(a)) | a G ^} C M''+^ . 

We assume that Pa is full-dimensional in that Mf^~^^ is its affine span. 

The upper facets of Pa are those facets whose outward pointing normal vector has positive 
last coordinate. Any face of an upper facet is an upper face. We illustrate this below when 
d = 1, where the displayed arrows are outward pointing normal vectors to upper facets. 




Figure 7. Upper faces and a regular subdivision. 



Projecting these upper facets to M*^ yields a regular polyhedral subdivision of the convex 
hull of A, which is the image of Pa. For our purposes, we will need to assume that the 
lifting function is generic in that all upper facets of Pa are simplices. In this case, we obtain 
a regular triangulation of A. This consists of a collection 

{Ai\i = l,...,m} 

of subsets of A, where each subset Ai consists of d+1 elements and spans a d-dimensional 
simplex. We regard all subsets of the facets Ai as faces of the triangulation. These simplices 
form a subdivision in that they cover the convex hull of A and any two with a non-empty 
intersection meet along a common face. 

The subdivision induced by the lifting function in Fig. [7] consists of three intervals result- 
ing from removal of the middle point of A, which does not participate in the subdivision 
as it is not lifted high enough. 

A set may have many regular triangulations, and not every point needs to participate in 
a given triangulation. Figure [8] shows the edges in four regular triangulations of 3 ^ fl Z^. 

Not every triangulation is regular. We may assume that a lifting function A for the 
triangulation of 4 ^ fl in Fig. [9] takes a constant value at the three interior points. The 
clockwise neighbor of any vertex of the big triangle must be lifted lower than that vertex. 
(Consider the figure they form with the parallel edge of the interior triangle.) Since the 
edge of the big triangle is lifted to a convex path, this is impossible, except in some M.C. 
Escher woodcuts. 

Definition 4.5. Let B = {ba | a G A} C be a collection of control points indexed by 
a finite set of exponents ^ C M'^ with d < n. Given a regular triangulation T = {Ai \ i = 
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Figure 8. Some regular triangulations. 




Figure 9. An irregular triangulation. 

1, . . . , m} of A we define the control polytope as follows. For each ci-simplex Ai in T, the 
corresponding points of B span a (possibly degenerate) simplex 

convjba I a G . 

The union of these simplices in R" forms the regular control polytope B{T) that is induced 
by the regular triangulation T of A. This is a simplicial complex in M" with vertices in B 
that has the same combinatorial type as the triangulation T of A. 

If the coordinate points (ea | a G ^) of are our control points (these are the vertices 
of ^■^), then the regular control polytope is just the geometric realization |T| of the 
simplicial complex T, which is a sub complex of the simplex /Sl'^. In general, B{T) is the 
image of this geometric realization |T| C under the projection vrg. 

Example 4.6. Let ^ := 3^ fl 1?, the exponents for a cubic Bezier triangle. Figure [TO] 
shows the three control polytopes corresponding to the last three regular triangulations of 
Figure [8|, all with the same control points. 

The reason that we introduce regular control polytopes is that they may be approximated 
by toric Bezier patches. 

Theorem 4.7. Let A C W^, w G , and B d he exponents, weights, and control 
points for a toric Bezier patch. Suppose that T is a regular triangulation of A induced by 
a lifting function X: A ^ M.. For each t > 0, let Ft: A ^ be the toric Bezier patch of 
shape A with control points B and weights t^^^^w^- Then, for any e > there exists a to 
such that ift > t^, the image Ft{A) lies within e of the control polytope B{T). 
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Figure 10. Three control polytopes. 



We prove Theorem 14.71 in Appendix [Bl Figure [TT] illustrates Theorem 14.71 for a cubic 
Bezier triangle with the control points of Example 14.61 The patch on the left is the cubic 
Bezier triangle with the weights of Example II. 5[ The second and third patches are its 
deformations corresponding to the lifting function inducing the leftmost control polytope 
of Fig. [TOl The values of t are 1,5, and 100, as we move from left to right. 




Figure 1 1 . Degeneration to the control polytope. 



An absolutely unpractical consequence of Theorem 14.71 is a universality result: Any 
surface which admits a triangulation that forms a regular control polytope may be approx- 
imated by a single Bezier patch. 

As with Theorem \4A\ the main idea behind the proof of Theorem 14.71 (which is given in 
Appendix [B]) is that for t large enough, the translated patch t.X^ C /^"^ can be made 
arbitrarily close to the geometric realization |T| C of the regular triangulation T. The 
result follows by projecting this into M" using tt^. 

In Appendix [B] we also prove a weak converse to Theorem 14. 7[ Namely if w.Xj\^ is 
sufficiently close to the geometric realization |T| of a triangulation T, then T is in fact the 
regular triangulation of A induced by the lifting function A (a) = log{wa)- 

Theorem 4.8. Let A C M'^ be a finite set of exponents. Suppose that \T\ C is the 
geometric realization of a triangulation T of A and there is a weight w such that the 
distance between Xjs^^yj and |T| is less than l/2{d+ 1). Then T is the regular triangulation 
induced by the lifting function A(a) = \og{wa). 
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Appendix A. Proofs of Injectivity Theorems 

Theorem 13.21 (Craciun-Feinberg) The map Gk is injective for every k G M'" if and only 
if the determinant of the Jacobian matrix, 



Jac(G 



k 



dxj 



does not vanish for any x E M" and any k G ^"^ 



Proof. We show the equivalence of the two statements, transforming one into the other. 
First, suppose there is a A; G M> so that Gk is not injective. Then there exist c, d G M" so 
that Gk{c) = Gk{d). Then we have 

m m 

= J2h{cy^-dy^)z, = ^fc^'((£)y>-i)z, , 

i=l i=l 

where k' G M> is defined by k'^ = kid^\ and | G M" has zth coordinate j-. In particular, 
Gk'i^) = Gk'iL), where i := (1, . . . , 1). Define v G M'^ by 

Vi := log(ci) - log(dj) , 
so that e^^'^ = {^Y', where ■ v is the Euclidean dot product, and we now have 

m 

(Al) = Y.K{ey'-~i)z, . 

1=1 

Define the univariate function by /(t) = (e* — l)/t for t 7^ and set /(O) = e. Then / is 
an increasing continuous bijection between M and M>. Define k" G M> by k" := k[f{yi ■ v). 
Then k'^ie^^'^ - 1) = Ki^i ' v), and fOTTD becomes 

m 

(A.2) = ^A;:'(y,-v)z, . 

i=l 

We claim that v lies in the kernel of the Jacobian matrix of Gk" evaluated at the point 
L. Indeed, let e^ := (0, . . . , 1, . . . , 0) G M" be the unit vector in the jth direction. Then 

where y^ = . . . , yi^n)- Since l^'-"^ = 1 and J2j yiVj = y,/ ■ v, we see that 

n n m m 

j=l ^ j=l i=l 1=1 

SO that V lies in the kernel of the Jacobian matrix of Gk" evaluated at l, which implies that 
the Jacobian determinant of Gk" vanishes at l. 

The theorem follows as these arguments are reversible. '^^^ 

Corollary 13.31 The map Gk (13.11) is injective for all k G M> if and only if [Yj ■ Zj) ■ [Yj ■ 

Zj) > for every I, J E ('™^) and at least one product Yj ■ Zj is non-zero. 
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Proof. Observe first that the Jacobian matrix Jac((j'fc) factors as the product of matrices 
S~^YDZ^, where 6 is the diagonal matrix with entries {xi, . . . ,Xn), D is the diagonal 
matrix with entries {kix^^ , . . . , kmX^"^) and Y and Z are the matrices whose columns 
are the vectors and Zj, respectively. If we apply the Binet-Cauchy Theorem to this 
factorization, we see that 

(A.3) detJac(G'fc) = x" ■ ^ W_kiX^' -Yi ■ Zi , 



/g/'MA i&I 



where i = fl, . . . , 11 



Suppose that (1/ ■ Zi) ■ {Yj ■ Zj) > for every I,J e C'^l) , and at least one product 
Yj ■ Zj is non-zero. Then all terms in the sum ( 1A.3I) have the same sign and not all are zero, 
and so the Jacobian does not vanish for any x G M" and k G M™. Thus Gk is injective for 
all G M>, by Theorem O 

Suppose that there are two subsets I, J G ('™') such that YjZj > and YjZj < 0. For 
teR> andK e ('^^l), define k{K,t) G by 



k{K,t), 



t ifjeK 
1 otherwise 



If we fix X G M.y, then the expansion (lA.Sp . implies that det Jac{Gk{K,t)){x) has the same 
sign as Y^Z^ when t S> 0, at least when Y^Z^ 7^ 0. 

We conclude that there is some k, x such that det Jac{Gk){x) > and some k, x such that 
det Jac(G'A;)(a;) < 0, and therefore some k,x such that det Jac(G'A;)(a;) = 0. This implies 
the corollary. '^^^ 



Appendix B. Three Toric Theorems 

Theorem 12.51 Suppose that A G M.'^ is a finite set of points with convex hull A. Let (3 = 
{/3a I a G A} be a collection of irrational toric Bezier functions for A. Then /3(A) = 
the closure of the image of ipj^. 



Proof. Let A° be the interior of A, which we assume has £ facets and is given by the facet 
inequalities < hi{x),i = 1, . . . ,i. Define two maps H : A° ^ and 4' '■ ~^ by 

H : X I — > {hi{x), hi{x)) , 

Then the map /? : A° ^ (whose image is dense in Xj) is the composition of the maps 
H and ip. 

Let us recall the definiton of the map (p^: iS^-^, 

LP A : ixi,...,Xd) I — > [x" : aeA] . 
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The theorem follows once we show that the map factors through the map ip_^. For this, 
define a new map /a : K> — > K> by 

e 

/a : u I — > t = (ti, . . . ,td) where tj := J^m^"" • 

1=1 

Then we claim that 

To see this, we compute the component of iplu) corresponding to a G 

i=l i=l i=l i=l 

Thus '?/'(«) = m'^V^^(/a(w)), as maps to M^. The common factor u'^ does not affect the 
image in iS'-^, which shows that ip = '^A ° /a and proves the theorem. ■'^^ 

A consequence of this proof of Theorem 12.51 is the derivation of equations which define 
the points of X^. This derivation is similar to, but easier than, the development of toric 
ideals in [191 Ch. 4], as we have monomials with arbitrary real-number exponents. 

Suppose that we have a linear relation among the points of A, 

(B.l) ^li^-a. = ^z^a ■ a , 

for some E M"^. Then the analytic binomial 



(B 2) n - n 



2^ =: 



aeA aeA 



vanishes on ipj^Mt,), considered as a point in M^. This follows from the easy calculation 

aeA 



Even after clearing denominators, the common zero set of the binomials (IB. 21) is not 
exactly the image ipj^^(M.y) in the simplex /S/-^, as the point fA{x) ^ is where the ray 

M> ■ (x^ I a G ^) C 

meets the simplex If we require that the binomial flB.2p is homogeneous in that 

^aG^ '"a ~ ^aGyt ^^cn it vauishes at every point of this ray and therefore on the image 
of V9_4 in zSZ"^. Since the coordinates are positive numbers, we may further assume that fIB.ip 
is an affine relation in that 

$^^a = 5^Z/a = 1 . 

aeA aeA 

These necessary conditions are also sufficient. 

Proposition B.3. A point z in ZS'-^ lies in Xjs^ if and only if we have 
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for all (j,,u eR-^ with Y^aeA 1^^ = 



aeA 



z/a = 1 and 



a 



• a 



ae^ aeA 



One way to see the sufficiency is to pick an affinely independent subset C oi A that 
affinely spans and use the formula x'^ = for a G C to solve for x G W^. Then the 
point z G satisfies this collection of binomials if and only ii z = This is also 

evident if we take logarithms of the coordinates. 

These arguments only work for points z in the interior of /S7-^. For points of X_4 with some 
coordinates zero, we use the recursive nature of polytopes and the toric Bezier functions. 
Namely, if we restrict the collection {(3^ \ a G A\ of toric Bezier functions to a face T 
of the convex hull of A, then those whose index a does not lie in T vanish, while those 
indexed by points of A lying in T specialize to toric Bezier functions for T. 

Theorem 14.71 Let A C , w G ]R>, and B C M" be exponents, weights, and control 
points for a toric Bezier patch. Suppose that T is a regular triangulation of A induced by 
a lifting function A: ^ — > R. For each t > 0, let Ft: A ^ be the toric Bezier patch of 
shape A with control points B and weights t^^^^w^- Then, for any e > there exists a to 
such that ift > to, the image Ft{A) lies within e of the control polytope B{T). 

Proof. The lifting function A : ^ ^ R inducing the triangulation T also induces an action of 
R> on where t G M> acts on a point z G by scaling its coordinates, {t.z)^, = t^^^^Zs,. 
Then Ft{A) is the image of t.X_4 under the projection vr^ : R". It suffices to show 

that t.X_A^w can be made arbtrarily close to \T\ C if we choose t large enough. 

We single out some equations from Proposition IB.3[ Suppose that a, b G ^ are points 
that do not lie in a common simplex of T. That is, the segment a, b is not an edge 
in the triangulation T, and therefore it meets the interior of some face of T so that 
there is a point common to the interiors of J-' and of a, b. (If a = b, so that a does not 
participate in the triangulation T, then this point is just a.) This gives the equality of 
convex combinations 



(B.4) 




where all coefficients are positive and J2c 



1. Thus 



holds on X^. The corresponding equation on is 



(B.5) Z^zl-'' = t/^a+(l-M)b-E.g^^c-c . "^a'^b . -Q ^u. 



Since a, b is not in the triangulation, points in the interior of the lifted segment 



(a,A(a)), (b,A(b)) 



lie below points of upper faces of the polytope Pa- We apply this observation to the 
point flB.4p . Its height in the lifted segment is /ia + (1— /i)b, while its height in the lift of 
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the face T is Xlcej^- t'c ■ c, and so 

/ia + (1— /i)b < 



This imphes that the exponent of t in fIB.SP is negative. Since the other terms on the right 
hand side are bounded, we see that the left hand side, and in fact the simple product Za,z\,^ 
may be made as small as we please by requiring that t be sufficiently large. 

Suppose that A consists of £ + (i + 1 elements. Repeating the previous argument for the 
(finitely many) pairs of points a, b which are not both in any simplex of T, we see that for 
any e > 0, there is a to such that if t > to and z G X_4 „,, then 

whenever a, b do not lie in a common simplex of T. In particular, at most one of Za or Zh 
can exceed e/2i. 

Let z G X_4 „,. Then there is some facet JF of T such that if a ^ JF, then < z^ < €/2i. 
Suppose that JF = {ao, ai, . . . , a^^} and set 



which is a point of the facet of the geometric realization |T| C Then 



\z - Zr\ 



Za^s 



^ai 



Za^e^ 



ai 



ZuLd^ad 



< 



aeA 

Za 



+ (1 



^ao 



< 2i 



2e 



as 1 = Ea ^a- 



We say that two subsets X and Y of Euclidean space are within a distance e if for every 
point X oi X there is some point y oiY within a distance e of x, and vice-versa. 

Theorem 14.81 Let A G M.'^ be a finite set of exponents. Suppose that \T\ C is 
the geometric realization of a triangulation T of A and there is a weight w such that the 
distance between X_4„, and \T\ is less than l/2{d+ 1). Then T is the regular triangulation 
induced by the lifting function A(a) = log(u;a). 

Proof. To show that T is the regular triangulation induced by the lifting function A whose 
value at a is log(wa), we must show that if a segment a, b between two points of A does 
not lie in the triangulation T, then its lift by A lies below the lift of some face JF of T. 
Set e := l/2{d+l). For each face JF of T, let xjr e be the barycenter of JF, 



E- 



aer 



where j^T is the number of points of A in JF, which is at most rf+l. If 2; is a point of Xj^^^ui 
within a distance e of Xjf, so that ||z — < e, then in particular no component of the 
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vector z — xjr has absolute value exceeding e. Thus we have the dichotomy 

< e = l/2(rf+l) if a^^, 

^ > z^ > l/#^-e > l/2(d+l) ifae^. 



Now suppose that the segment a, b does not lie in the triangulation T. Then there is 
a face of the triangulation whose interior meets the interior of this segment. That is, 
there is an equality of convex combinations (lB.4p and a corresponding equation that holds 
for points z G X_4^, 

4^1' ■ n < = ■ n <^ • 

Suppose that 2; is a point of X_4 „, that lies within a distance of l/2{d+l) of the barycenter 
of the face J-'. Then, by the estimates flB.6p . we have 

TT > z^z}-" ■ TT Wc'^ = Kw}-'' ■ n <^ > w^w}-" ■ ,] , . 

2 J-1 ^ a b 11 c a b 11 c a b 2 

canceling the common factor of l/2((i+l) and taking logarithms, we obtain 

log(Wc) > /ilog(Wa) + (1 - yU) log(t(7b) , 



which implies that the point ( IB. 41) common to the segment a, b and the face JF of T is 



lifted higher in the face T than in the segment a, b, and so the lift of the segment a, b by 
A lies below the lift of the face JF. As this is true for all segments, we see that T is the 
triangulation induced by the lifting function A. 
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